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Abstract 

Using the QCD hybrid (moments-Laplace) sum rule, we show semi- analytically 
that, in the limit M& — > oo, the q 2 and M& behaviours of the heavy-to-light 
exclusive (B — > p (n) semileptonic as well as the B — > /ry rare) decay- 
form factors are universally dominated by the contribution of the soft light- 
quark condensate rather than that of the hard perturbative diagram. The 
QCD-analytic q 2 behaviour of the form factors is a polynomial in q 2 /M£, 
which mimics quite well the usual pole parametrization, except in the case 
of the Af form factor, where there is a significant deviation from this po- 
lar form. The M^-dependence of the form factors expected from HQET 
and lattice results is recovered. We extract with a good accuracy the ratios: 
V B (0)/Af (0) ~ Af (0)/A? (0) ~ 1.11±0.01, and Af(0)/Ff(0) ~ 1.18±0.06; 



combined with the "world average" value of f? (0) or/and Ff(0), these ra- 
tios lead to the decay rates: Tg^ weP ~ (4.3 ± 0.7) x \V ub \ 2 x 10 12 s _1 , 
^B^peu/^B~*neu — -9 ± 0.2, and to the ratios of the p— polarised rates: 
r+/r_ ~ 0.20 ± 0.01, a = 2V L /T T - 1 ~ -(0.60 ± 0.01). 
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1 Introduction 



In previous papers |I[ |2|, we have introduced the hybrid (moments-Laplace) 
sum rule (HSR), which is more appropriate than the popular double expo- 
nential Laplace (Borel) sum rule (DLSR) for studying the form factors of 
a heavy-to-light quark transition; indeed, the hybrid sum rule has a well- 
defined behaviour when the heavy quark mass tends to infinity. In |||, we 
studied analytically with the HSR the M^-dependence of the B — ► K*^f form 
factor and found that it is dominated by the light-quark condensate and be- 
haves like \[Mb at q 2 = 0. We have also noticed in |j| that the light-quark 
condensate effect is important in the numerical evaluation of the B — > p (71) 
semileptonic form factors, while it has been noticed numerically in || using 
the DLSR that for the B — ^ p semi-leptonic decays, the q 2 behaviour of 
the Af form factor in the time-like region is very different from the one ex- 
pected from the standard pole representation. In this paper, we shall study 
analytically the M^-behaviour of the different form factors for a better under- 
standing of the previous numerical observations. As a consequence, we shall 
re-examine with our analytic expression the validity of the g 2 -dependence ob- 
tained numerically in [[3], although we shall mainly concentrate our analysis 
in the Euclidian region (q 2 < 0). There, the QCD calculations of the three- 
point function are reliable; also the lattice results have more data points. For 
this purpose, we shall analyse the form factors of the B — > n(p) semileptonic 
and B — > /ry rare processes defined in a standard way as: 



(p(p')\u lfM (l- l5 )b\B(p)) = (M B + M p )A^ - ^ — e*p'(p + p' „ 

2V 



M B + M P 

(n(p')\u lfl b\B(p)) = / + (p+p% + /-(P-PV 
<p(p')\sa, u (^^q v b\B(p)> = *We>V ff ^ P 

+ {^M 2 B -M 2 )-e*q(p + p'),} 
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In the QCD spectral sum rules (QSSR) evaluation of the form factors, we 
shall consider the generic three-point function: 

V{p,p',q) = -Jd 4 xJ d A ywp{ip'x - ipy)(0\TJ L {x)O{0)J b {y)\0), (2) 

whose Lorentz decompositions are analogous to the previous hadronic ampli- 
tudes. Here Jl = ujpd (J L = {m u + rridjui^^d) is the bilinear quark current 
having the quantum numbers of the p (tt) mesons; J b = (M b + m d )di^^b is 
the quark current associated to the S-meson; O = try^u is the charged weak 
current for the semileptonic transition, while O = b^a^ v q v is the penguin 
operator for the rare decay. The vertex function obeys the double dispersion 
relation: 

00 ds r°° ds' 



V(p 2 ,p' 2 ,q 2 ) = — -J -ImV(s,s',q 2 ) + ... (3) 

Air z Jm? s — p Jo s — p' z 
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As already emphasized in @, we shall work with the HSR: 



d P 2 J P 2 =0 v v v ' 

1 /"DO f7e /"OO 

?L7^Io ds'eM-r's')ImV(s,s', q % (4) 



rather than with the DLSR (C is the Laplace transform operator). This sum 
rule guarantees that terms of the type: 

Ml 1 



— , (5) 
(M 2 -]9 2 )V 2fe 

which appear in the successive evaluation of the Wilson coefficients of high- 
dimension operators, will not spoil the OPE for M b —>■ oo unlike the case of 
the double Laplace transform sum rule, which blows up in this limit for some 
of its applications in the heavy-to-light transitions. 

In order to come to observables, we insert intermediate states between the 
charged weak and hadronic currents in (2), while we smear the higher-states 
effects with the discontinuity of the QCD graphs from a threshold t c (t' c ) for 
the heavy (light) mesons. Therefore, we have the sum rule: 

Hv\j * 2C L f B ^±exp(-Mlr) 
1 r*c ds K 



4tt 2 



rt c rig rt c 

/ 9 ^tt/ ds'expi-Ts'^mVpT^s'^q^+NPT. (6) 

J Ml S n+L JO 

PT (NPT) refers to perturbative (non-perturbative) contributions; Cl = 
fpMp for light pseudoscalar mesons, while Cl = My/(2jv) for light vector 
mesons; Mi is the light meson mass. The decay constants are normalized as: 

K + M Q ) (0\q{i l5 )Q\P) = V2M 2 P f P 

(0\qj„Q\V) = e;V2^. (7) 

F(q 2 ) is the form factor of interest. For our purpose, we shall consider the 
expression of the decay constant fg from moments sum rule at the same 
order (i.e. to leading order) [|J: 

2/1 _ 3 2 ftc ds (s-M 2 ) 2 <qq> f n 2 (n 2 + 1) ( M 2 \] 



(M 2 B ) n2 ~ l 8tt 2 b J Ml s" 2+1 s M% n2 ~ l \ 4 V M 6/J 

(8) 

For convenience, we shall work with the non-relativistic energy parameters 
E and 5M {b y. 

s = (M b + Ef and SM (b) = M B - M b , (9) 

where, as we saw in the analysis of the two-point correlator, the continuum 
energy E c is [|J: 

E? ~ (1.08 ± 0.26) GeV 
#f ~ (1.30 ± 0.10) GeV 

£ c °° ~ (1.5 ~ 1.7) GeV. (10) 
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In terms of these continuum energies, and at large values of M b , the decay 
constant reads J3J: 



+- (2n 2 + 3)(n 2 + l) + - -f- 



'3fV 7T 2 (<?<?> A n 2 (n 2 + l)M 2 



M 6 ; 2 (Efy v 4 M 



2 The _B — ► p semileptonic decay 

The corresponding form factors defined in (1) have been estimated with the 
HSR Jl| and the DLSR 0, 0. Instead of taking the average values from the 
two methods as was done in [J], we shall only consider the HSR estimates, 
because of the drawbacks previously found in the DLSR approach: 

Af(0) ~ 0.16 -0.41, A%(0) ~ 0.26- 0.58, V B (0) ~ 0.28 - 0.61. (12) 

The errors in these numbers are large, as the HSR has no n-stability. In the 
following, we derive semi-analytic formulae for the form factors. Using the 
leading order in a s QCD results of the three-point function, and including 
the effect of the dimension-5 operators as given in || , one deduces the sum 
rule (q 2 < 0): 



<^4<<(fH i -&+* ,5,+ ^ (i3) 



with: 



_ ( lp\ M b 2 



5 (5) 



r'M 2 j 1^3 3 2 
' n K 1 n n 



6 I " t'MI \~ 4" 4' 



g ^(n + 1) ( ~n - 1 ) + 2r'M b 2 (l + 2n) + 2(n + l)gV 



M? V y V2 



(14) 



where is the integral from the perturbative expression of the spectral 
function. It is constant for Mb — > oo. Its value and behaviour at finite values 
of Mb and for q 2 = is given in Fig. 1. At M& = 4.6 GeV, it reads: ~ 
(9./t oo.e) GeV e and behaves to leading order in 1/M& as t' 2 E c / (qq), which 
is reassuring as it gives a clear meaning of the expansion in (13). The other 
values of the QCD parameters are §: (qq) = -(189 MeV) 3 (log M b /A) 12/23 
and M 2 = (0.80 ± 0.01) GeV 2 from the analysis of the B, B* sum rules. The 
p-meson coupling is 7 P ~ 2.55. 

One can deduce from the previous expression that Af is dominated by the 
light-quark condensate in the l/M&-expansion counting rule. Moreover, the 
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perturbative contribution is also numerically small at the 6-mass. The ab- 
sence of the n-stability is explicit from our formula, due to the meson-quark 
mass difference entering the overall factor. This effect could be however 
minimized by using the expression of Jb in (11) and by imposing that the 
effects due to the meson-quark mass differences from the three- and two- 
point functions compensate each other to leading order. This is realized by 
choosing: 

2n = n 2 -^, (15) 

which, fixes n to be about 2, in view of the fact that the two-point function 
stabilizes for n 2 ~4-5. In this way, one would obtain the leading-order result 
in a s \ 

Af~0.3-0.6, (16) 

where we have used the leading-order value fg° ~ 1.24/^. However, al- 
though this result is consistent with previous numerical fits in (12) and in 
@, we only consider it as an indication of a consistency rather than a safe 
estimate because of the previous drawbacks for the n— stability. One should 
also keep it in mind that the values given in (12) correspond to the value 
of fs — 1-6/tt, which includes the radiative corrections of the two-point cor- 
relator and which corresponds to smaller values of n. Improvements of the 
result in (16) need (of course) an evaluation of the radiative corrections for 
the three-point function. The g 2 -dependence of Af can be obtained with 
good accuracy, without imposing the previous constraint. We obtain the nu- 
merical result in Fig. 2, which is well approximated by the effect from the 
light-quark condensate alone: 

T?B( a 2\ = A l^) „ I _ ( 17 ) 

Performing an analytic continuation of this result in the time-like region, we 
reproduce the numerical result from the DLSR @(see Fig. 2), which indicates 
that the result is independent of the form of the sum rule used, while in the 
time-like region the perturbative contribution still remains a small correction 
of the light-quark condensate one. This result is clearly in contradiction with 
the standard pole-dominance parametrization, as, indeed, the form factor 
decreases for increasing q 2 - values. A test of this result needs improved lattice 
measurements over the ones available in 0. From the previous expressions, 
and using the fact that fs behaves as 1/y/Mf,, one can also predict the M b - 
behaviour of the form factor at q^ax — + 2M p M;,: 



2 N 

1 r\j 



maxJ nr, 



in accordance with the expectations from HQET || and the lattice results 
0. The analysis of the V B and Af form factors will be done in the same 
way. Here, one can realize that the inclusion of the higher dimension-5 and -6 
condensates tends to destabilize the results, although these still remain small 
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corrections to the leading-order results. Then, neglecting these destabilizing 
terms, one has: 




V B (q 2 ) = 

Ai(q 2 ) * (W) ^ (^) <! 1 I -?*+...} (19) 



2 /. 



with: 




M 



ft = Its] «^e,p,M>'). (20) 

Zy g are integrals from the perturbative spectral functions, which also be- 
have like Tqo to leading order in 1/M&. They are given in Fig. 1 for 
q 2 = and for different values of M&. As expected, they are constant when 
M b — > 00, although, as in the previous case, the asymptotic limit is reached 
very slowly. Here, the n-stability of the analysis is also destroyed by the 
overall (Mb/ M\,) 2n factor, which hopefully disappears when we work with 
the ratios of form factors. We show in Fig. 2 the g 2 -dependence of the 
normalized V B and A B , which is very weak since the dominant light-quark 
condensate contribution has no g 2 -dependence. The small increase with q 2 is 
due to the g 2 -dependence of the small and non-leading contribution from the 
perturbative graph. Lattice points in the Euclidian g 2 -region [7|] agree with 
our results. An analytic continuation of our results at time-like q 2 agrees 
qualitatively with the one in ||. The numerical difference in this region 
is due to the relative increase of the perturbative contribution in the time- 
like region due to the effect of the additional non-Landau-type singularities. 
However, this effect does not influence the M b behaviour of the form factors 
at q^axy which can be safely obtained from the leading-order expression given 
by the light-quark condensate. One can deduce: 

,2 \ . [77 aB,2 



V B (qL x ) ~ V M fc> A 2 (iLx) ~ \JM b . (21) 

This result is in agreement with HQET and lattice data points. Finally, we 
can also extract the ratios of form factors. At the r'-maxima and at the 
n-maxima or inflexion point, we obtain from Fig. 3: 

r ^iM" r ^^oy" 1 - 11±0 - 01 ' (22) 

where the accuracy is obviously due to the cancellation of systematics in the 
ratios. This result is again consistent with the lattice results J7|, but more 
accurate. 
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3 The B — > 7r semileptonic decay 

The relevant form factor defined in (1) has been numerically estimated within 
the HSR with the value 0: 

/f (0) ~ 0.20 ± 0.05, (23) 

where the contribution of the 7r'(1.3) meson has been included for improving 
the sum rule variable stability of the result. In this paper, we propose to 
explain the meaning of this numerical result from an analytic expression 
of the sum rule. Using the QCD expression given in 0, we obtain for a 
pseudoscalar current describing the pion: 

2n 



4/.m2 f B \M b J { Mi 



(24) 



where X n is the spectral integral coming from the perturbative graph. Its 
value at q 2 = for different values of M& is shown in Fig. 1. It indicates 
that at M b = 4.6 GeV, the perturbative contribution, although large, still 
remains a correction compared with the light-quark condensate term; 6^ is 
the correction due to the dimension-5 condensate and reads: 



^-^HS^ 1 ^"- 1 )}' (25) 

One can use the well-known PCAC relation 

(m u + m d ) (qq) = -m 2 J 2 , f v = 93.3 MeV (26) 

into the previous sum rule in order to express f? in terms of the meson 
couplings. Unlike the case of the B — > p form factors where the scale depen- 
dence is contained in (qq), /+ is manifestly renormalization-group-invariant. 
It should be noted, as in the case of the sum rule determination of the cupTT 
coupling , that the /^-dependence appears indirectly via (26) in a correla- 
tor evaluated in the deep Euclidian region, while the pion is off shell, which 
is quite different from soft-pion techniques with an on-shell Goldstone bo- 
son. One can also deduce from (24) that for large Mb, f? behaves like a/M&. 
In this limit the g 2 -dependence is rather weak, as it comes only from the 
non-leading 1/M& contributions; we therefore have, to a good accuracy: 

f+(q 2 max ) * /» ~ y/m- (27) 

As in the previous case, the slight difference between the g 2 -behaviour in 
the time-like region and the one from that obtained in ||, at a finite value 
of Mb(=4.6 GeV), is only due to a numerical enhancement caused by the 
non-Landau singularities of the perturbative contribution in this region, but 
does not disturb the M&-behaviour of the form factor. Finally, we extract 
the ratio of the form factor: 

r. - J^- (28) 

Unfortunately, we do not have stabilities, as the stability points are different 
for each form factor, which is mainly due to the huge mass- difference between 
the p and 7r mesons. 
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4 The B — ► p7 rare decay 

We can use the previous results into the HQET |Sj relation among the differ- 
ent form factors of the rare B — > p7 decay (Ff = F^ p ) and the semileptonic 
ones. This relation reads around q^ax '■ 

Fl (g } " 2M^ M B + M P + 2M B ^ (g } ' (29) 

from which we deduce: 

if (<?LJ ~ \M- (30) 

However, we can also study, directly from the sum rule, the g 2 -dependence of 
Ff. Using the fact that the corresponding sum rule is also dominated by the 
light-quark condensate for Mb — ► oo 0, an evaluation of this contribution, 
at q 2 7^ 0, shows that the light-quark condensate effect has no g 2 -dependence 
to leading order. Then, we can deduce, to a good accuracy: 

F l{<lmax) ^ *f (0) ~ \[M b . (31) 

Let us now come back to the parametrization of the form factor at q 2 = 
0. We have given in an expanded interpolating formula that involves 
1 /Mb and 1/M 2 corrections due to the meson-quark mass difference, to fs 
and to higher- dimension condensates. Here, we present a slightly modified 
expression, which is: 



2 / 



B 




FfWc-Uqq)^ (^) <1 + £L+ L (32) 



with: 



ft) exp(M ^'» 



p 7 

J 7 ~ (20 ± 4) GeV 2 for M b > 4.6 GeV, (33) 

where we have neglected the effects of higher-dimension condensates; Z 7 is 
the perturbative spectral integral. One should notice that unlike the other 
spectral integrals in Fig. 1, Z 7 reaches quickly the asymptotic limit when 
Mb — * oo. Using the estimated value of Ff(0) in [|J, we can have, in units 
of GeV: 

d 1.6 x 1(T 2 / 20±4\ 
F i = 7 1 + ^5- . (34) 



which leads of course to the same formula at large Mb as in ||. However, due 
to the large coefficient of the perturbative contribution, it indicates that an 
extrapolation of the result obtained at low values of M c is quite dangerous, 
as it may lead to a wrong M^-behaviour of the form factor at large mass. 
One should notice that (34) and the one in || lead to the same numerical 
value of F 1 Z3 (0). Proceeding as for the former cases, we can also extract the 
ratio: 

AM 
Ff(0) 

from the analysis of the r'- and n-stability shown in Fig. 3. 



r 7 = -=?±+ ~ 1.18 ±0.06, (35) 
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5 Values of the _B-form factors 



The safest prediction of the absolute value of the form factors available at 
present, where different versions of the sum rules and lattice calculations 
have a consensus, is the one for /+ (0): 

/£ (0) ~ 0.26 ± 0.12 ± 0.04 Lattice [7] 

0.26 ± 0.03 DLSR ||(see alsog) 

0.23 ±0.02 HSR+DLSR [HI 



0.27 ±0.03 Light-cone [10 



(36) 

from which one can deduce the "world average" : 

/£ (0) ~ 0.25 ±0.02. (37) 

For estimating Af(0), one can use the present most reliable estimate of Ff 
i, 0: 

F?(0) ~ 0.27 ±0.03, (38) 

where we have used the strength of the 5 , [/(3)-breakings obtained in [fj], in 
order to convert the result for B — > K*^ of (ll| into the B — > p7 of interest 
here. Then, we deduce: 

Af(0) ~ 0.32 ±0.02, (39) 

which is consistent with a direct estimate [[l], but the result is again more 
accurate. 



6 .B-semileptonic-decay rates 

We are now in a good position to predict the different decay rates. In so 
doing, we shall use the pole parametrization, except for the Af form factor. 
For the B — > 7r, we shall use the experimental value 5.32 GeV of the B* mass. 
For the B — > p, we shall use the fitted value (6.6 ± 0.6) GeV || for the pole 
mass associated to A® and V B . For Af , we use the linear form suggested 
by (13), with an effective mass of (5.3 ± 0.7) GeV, which we have adjusted 
from the numerical behaviour given in || (we have not tried to reproduce 
the change of the behaviours for t ~ (0.76 — 0.95)Af| obtained in 0, which 
is a minor effect). Using the standard definitions and notations, we obtain: 

~ (3.6 ± 0.6) x \V ub \ 2 x 10 12 s- 1 (40) 

We also obtain the following ratios: 

B^peu ^ Q 9±0 2; _t ~ o.20±0.01, a = 2-^-1- -(0.60 ±0.01). 

1 B-^weu 1 - 1 T 

(41) 
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Thanks to a better control of the ratios of form factors, the ratio of the B 
decays into n over the p can be predicted, to a good accuracy. It becomes 
compatible with the prediction obtained by only retaining the contribution 
of the vector component of the form factors. Our predictions are compatible 
with the ones in || except for r + /T_, where the one in || is about one order 
of magnitude smaller. The difference of two of these three quantities with 
ones in QXJ] (the large branching ratio into p over 7r and the positive value of 
the asymmetry a in and in most other pole dominance models for Af ) is 
mainly due to the different q 2 — behaviour of Af used here. 

7 Conclusions 

We have studied, using the QCD hybrid sum rule, the Mb- and g 2 -behaviours 
of the heavy-to-light transition form factors. We find that these quantities are 
dominated in a universal way by the light-quark condensate contribution. 
The Mfe-dependence obtained here is in perfect agreement with the HQET 
and lattice expectations. 

The g 2 -dependence of the Af form factor, which is mainly due to the one from 
the light-quark condensate contribution, is in clear contradiction with the one 
expected from a pole parametrization. The other form factors can mimic 
numerically this pole parametrization. Our QCD-analytic g 2 -behaviours 
confirm the previous numerical results given in ||. 

We have also shown that it can be incorrect to derive the M^-behaviour of 
the form factors at q 2 = by combining the HQET result at g^ aa . with the 
pole parametrization. 

We have also shown that the unusual q 2 — behaviour of the Af form fac- 
tor affects strongly the branching ratio of B — > p over B — ► 7r and the p- 
polarisation parameter a. A measurement of these quantities complemented 
by the one of the q 2 — behaviour of the form factor should provide a good test 
of the sum rules approach. 

We want also to stress that the extrapolation of the results obtained in this 
paper to the case of the D-meson would be too audacious: the uses of the 
HSR in that case cannot be rigorously justified since the value of the c-quark 
mass is smaller, although it may lead to acceptable phenomenological results. 
We are investigating this point at present. 
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Figure captions 

Fig. 1 Mb-dependence of the perturbative spectral integrals at q 2 = 0. 
Fig. 2 g 2 -behaviour of the normalized form factors: R\ = Af (q 2 )/Af (0), 
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R 2 ee A B (q 2 )/A B (0), R v = V B (q*)/V B (Q) and R w = f B (q 2 )/f B (0). The 
squared points in the timelike region are from 0. 

Fig. 3 t'- and n-dependences of the ratios of form factors at q 2 = 0: 
r 2 = A B (0)/Af (0), r v = V B (0)/Af (0) and r 7 = (0)/Ff (0). 
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